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Abstract 

The existing experimental evidence on the vector and axial-vector coupling 
constants, Gy and Ga-, for neutron d-decay exhibits two prominent problems; 
(i) the unitarity for the first row of the CKM matrix seems violated; (ii) one 
obtains different values of GaIGv according to whether one uses as input 
the neutron lifetime or neutron-decay correlation observables. We show that 
the in-medium modihcation of Ga can influence both the inner and outer 
radiative corrections used in deducing Gy from the super-allowed Fermi tran¬ 
sitions and that this effect may resolve the long-standing problems (i) and (ii) 

simultaneously. 
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The vector coupling constant Gy and the axial-vector coupling constant Ga for neutron 
/d-decay are quantities of primary importance that have been studied extensively. The 
recent remarkable progress in the determination of Gy and Ga highlights two very notable 
problems. One problem is concerned with Gy and the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element Vud- Comparison of Gy with obtained from the p-decay rate gives the 
empirical value of Vud, and the existing data give |I| 

Kd = 0.9740 ± 0.0005. (1) 

Since in the standard model the CKM matrix is unitary, the hrst row of the CKM matrix 
must satisfy |14dp -|- |14sP -1- |14feP = 1- However, Vud in Eq.(2) combined with the available 
information on Vus and Vub leads to 

iKdl" + \Vus\^ + iKfeP = 0.9972 ± 0.0013, (2) 

indicating a violation of unitarity at the 98% conhdence level [Q. Since this violation, 
if conhrmed, constitutes a clear signal for physics beyond the standard model, and since 
the contribution of Vud dominates the sum in Eq.(g), it is vitally important to assess the 
reliability of Vud in Eq. (|^) . 

The second problem concerns A = GA/Gy. The neutron-decay angular distribution 
experiments give P| 

Aa = 1.2601 ±0.0025, (3) 

whereas the neutron lifetime measurements, which have recently greatly improved, give 0 

A^ = 1.2681 ±0.0033. (4) 

(For the neutron lifetime we use the value 888 ±3 sec which is between the two recommended 
values 887 ± 2.0 sec and 889.2 ± 2.2 sec.) The origin of this discrepancy in A is at present 
unknown, a situation that is rather disturbing. 

In this Letter we point out the possibility that the above two problems are related to 
each other and that they may be solved simultaneously by taking into account certain types 
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of hadronic effects which have not been considered so far. As is well known, Gy is obtained 
from the observed ft values for O’*" ^ O’*" nuclear /d-decays after applying radiative corrections 
and isospin-mixing corrections. Although these corrections have been studied in great detail 
by many authors, we present here a new type of nuclear correction which increases Gy, and 
hence Vud as well, in such a way that the unitarity condition is satished. This increase in 
Gy decreases the value of A,- due to the fact that the neutron decay rate is proportional to 
Gy + ?iG\ = Gy(l-|-3A^). By contrast, A^ remains unaffected, because it is directly related 
to the correlation observables in neutron /3-decay and hence has no connection with the ft 
values of the O’*" —> O’*" transitions. We will argue that our new value of A,- is compatible 
with Aa- 

In deducing the empirical value of Vud from comparison of the rates of O’*" —O’*" nuclear 
/3-decays and the strength of the pure leptonic /x-decay, an accurate estimate of radiative 
corrections is the most important step (see for example, [|l],^|-^ and references therein). The 
bench mark work which constitutes the basis of all the recent analyses is due to Marciano 
and Sirlin |^, who hnd the following corrections 

(y 

^ [9{Em) + 41n (mz/M) + In (M/Ma) + 2G + A,]. (5) 

The hrst term, often called the outer correction, represents a spectrum-average effect which 
depends on the /3-ray end-point energy Em, and its contribution amounts to 0.01 i. 
Of the remaining terms (inner corrections), the 4In (m^/M) term represents the dominant 
model-independent short-distance contribution (M = proton mass, mz = Z-boson mass), 
whereas In (M/Ma) + 20 are axial-current induced contributions. M^(~l GeV) is a low- 
energy cutoff characterizing the short-distance part of the yhh box diagram, while G repre¬ 
sents the long-distance correction. Numerically, the contributions of the short-distance and 
long-distance parts are comparable, each contributing about 0.001 (see, e.g. ref. 0). Fi¬ 
nally, Ag is a perturbative QCD correction, which turns out to be very small. The standard 
radiative corrections above have been obtained for the weak quark current g 7 ^(l -|- 75 ) 5 . In 
applying theses results to nucleon and nuclear /3-decays, one needs to translate the quark- 


3 


based description to hadronic descriptions, and this step is considered to introduce a certain 
amount of model dependence. 

The fact that the axial current can contribute to the radiative correction of the Fermi 
transition means that, if one uses the nucleon weak current instead of the quark current, 
the results will in general depend on A, the strength of the axial-current coupling, and that 
the well-established nuclear-medium dependence of A is capable of affecting the radiative 
correction. In fact, parts of such effects have been studied in the literature. Several authors 
1^-^ investigated the two-nucleon contribution to the C term in Eq. (|^), while Towner P] 
estimated the influence of the in-medium modihcation of A on the C term. Changes in C 
due to these effects were reported to be small. The existing work is limited to the C term 
because in the treatment of Ref. |^, the ln(M/M^) term, representing the short-distance 
effects was considered to be independent of hadron/nuclear models, leaving only the C term 
as a mo del-dependent contribution. We note, however, that, although the separation of the 
axial-vector induced corrections into short- and long-distance parts has a clear-cut meaning 
in the quark picture, there is a possibility that this separation becomes less well dehned 
when one switches (as one must do at some point) to hadronic descriptions. 

In this note we estimate the effects of the in-medium modihcation of A, using the nucleonic 
weak current from the outset, and without invoking separation into short- and long-distance 
contributions. This means that our calculation in principle incorporates the combined effects 
of the ln(M/M^) and C terms, even though we do not try to make direct contact with 
the quark-model-based classihcation used in Eq. (j^). As will be demonstrated below, 
our approach indeed gives rise to a cut-off dependent term which is a function of A and 
hence can be inhuenced by nuclear-medium effects. It is important to emphasize, however, 
that we work here with the nucleonic weak current instead of the quark weak current not 
because we think the former has more fundamental meaning. Rather, the reason is that the 
nucleonic current description provides an interesting alternative framework for incorporating 
the nuclear medium effects. 

Our work has some overlap with that of Towner P] in that the effects of in-medium 
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modification of A are considered but, as mentioned above, the calculation in ||^ was limited 
to the C term. Obviously, a calculation of the entire A-dependent effect (rather than just 
the C term) is an extremely difficult task. In the present exploratory work, therefore, we 
present an estimate based on schematic models. Thus, we use the nucleon weak current 

+ A75)m„, (6) 

and assume that the nuclear-medium effect can be represented by changing the free space 
value A into an in-medium value A*. As for electromagnetic interactions we take into account 
the nucleon charge and magnetic moments and the electron charge. To calculate the loop 
diagrams, we employ two simplihed approaches. In the hrst approach we use a single common 
cut-off parameter Aq instead of form factors at the vertices. In the second approach we use 
a common form factor A^/(g^ -|- A^), both at the weak and at the electromagnetic vertices, 
where q is transferred momentum, and A is the parameter of the same order of magnitude 
as the nucleon mass M. 

For each of the two cases, a simple but quite lengthy calculation leads to the following 
expression for a A-dependent contribution to the inner radiative correction: 

Case I: 

A„ = a|^^ [iJIp - 7lJ(eo - 1) In (1 - eo) + 2[eo - 2 - In (1 - ^o)]) , (7) 

Case II; 

= A^ |^(7Ip-7 lJ(^-^lnO + ^ + 2^1 n^+^ + ^ln^-^ln^^ , (8) 

where /Ip and /!„ are the proton and neutron anomalous magnetic moments, = M^/Aq and 
^ = M^/A^. We have checked that the two rather different phenomenological approaches 
used here give numerically similar results within reasonable variations of the model param¬ 
eters, Aq and A. Even in the absence of these cut-off parameters (point-like nucleons with 
a cut-off at the scale of the vector boson mass) the calculated inner corrections does not 
exceed the values in Eqs.(0) and (|) by more than an order of magnitude. This indicates a 
certain degree of stability of our results. 
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In the above we concentrated on the inner corrections. We now turn our attention to 
the outer correction, which is dehned as an ii^m-dependent correction. As is well known [Q, 
if one only retains terms of 0{a), ignoring terms of 0{aE/M) and 0{a{E/M)ln{M/E)) 
{E= electron energy), then the outer correction becomes model-independent and hence also 
A-independent. The term 'g[Em) in Eq. represents this model-independent contribu¬ 
tion. However, at the level of precision we are interested in, it becomes necessary to retain 
terms of 0{aE/M) and 0{a{E/M) \n{M/E)), and these terms lead to a A-dependent outer 
correction. We write a generic form of the A-dependent outer correction as 

Pouter = (9) 

The explicit forms of /(E^), which depend on the models used, are very lengthy and will be 
given elsewhere. We remark that the separation of the radiative corrections into the terms 
in Eq.(^ is specihc to the quark model, and therefore it is not possible to establish termwise 
correspondence between the terms in Eq.(H) and those in Eqs.(|^ - (^. 

We denote by the leading-order A-dependent radiative correction, which is the sum 
of the inner correction in Eq. (|^ or Eq. (^), and the outer correction in Eq. (^: 

(y 

= dinner ^outer — ^ 

Ztt 

Here J contains information about the nucleon structure including its weak and electromag¬ 
netic form factors. Assuming that the cut-off parameters in Eqs.(0)-(^) are controlled by a 
strong interactions scale (Aq ~ A ~ M) and that a characteristic energy of nuclear /3-decay 
is Eyyi 5 MeV, the value of J turns out to be typically 5 ~ 25. In the subsequent numerics, 
we use J = 15 ± 10. 

Since we have a general expression (within our schematic model) for the A-dependent part 
of the 0{a) radiative corrections, we may use it to study the influences of the in-medium 
modihcation of A. We expect that, if a description based on the nucleon weak current 
achieves a sufficient level of precision, then the result corresponding to the case with A fixed 
at the free nucleon value should essentially reproduce the results of the existing analyses 
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of the relevant observables based on Eq.(^. Then we may assnme that small pertnrbative 
changes in the A-dependent radiative correction due to in-medium modihcation of A can be 
estimated with the present schematic model. We use Eq.([T0|) in this limited context. 

We introduce the in-medium renormalization parameter p as p = A*/A. For a wide range 
of the periodic table, a typical value of p is p ~ 0. 8 |igg. For our present purpose, we £x 
p at 0.8. Then, by taking the difference of Aa, Eq. ( 0 ), between the nuclear and free cases, 
we obtain 

A = Aa-Aa. = A(l-p);^J. (11) 

ZTT 

This leads to a corrected value of the CKM matrix element Vud 


Vud Wd(l + A/2), 


( 12 ) 


and, since the probability of neutron /3-decay is proportional to the |14dP, to a corrected 
expression for the neutron life time 


= 


2{Ft) (1-A) 


(13) 


/Kln2(l + 3A2) 

Here F and /r are nuclear and neutron phase space factors [§|,0. Eqs. (P!2|) and (P!B|) sum¬ 
marize how the medium dependence of A encoded in A can change the deduction of Vud 
and the analysis of r„. For the choice of the values of parameters explained above, we have 
A = 0.0044 ±0.0029. 


The new formula Eq.(^) changes the value of Vud from Vud = 0.9740 ± 0.0005 to Vud = 
0.9761 ± 0.0015. The error for the CKM matrix element has been estimated simply as the 
sum of the squared errors, one coming from Eq.(|ID and the other coming from the allowed 
range of J. The new value of the CKM matrix element in turn leads to 


\Vud? ± IK.P ± \Vub? = 1.0013 ± 0.0031 


(14) 


in good agreement with the unitarity condition. Moreover, the use of Eq. (pAf ) with A = 
0.0044 ± 0.0029 decreases A,- from the value in Eq.(|p to At- = 1.2655 ± 0.0034. This new 
value agrees with A^i, Eq.(||) within the error bars. 


7 



Thus, we have illustrated the possibility that the inclusion of the effect of the in-medium 
modihcation of in the entire radiative correction (not only in the C term) is capable of 
removing the two outstanding problems regarding the basic coupling constants Gy and Ga- 
Our present treatment is admittedly very schematic and subject to various improvements, 
and the particular numerical successes reported above are simply indication of the possible 
effect. However, we believe that our calculation does indicates the importance of the new 
type of radiative correction considered here and that future analyses of the O"*" —O’*" /3- 
transitions need to include the effects discussed in this article. 
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